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Abstract: We discuss the convergence of level truncation in bosonic open string field 



QJ , theory. As a test case we consider the calculation of the quartic tachyon coupling 74. We 

' determine the exact contribution from states up to level L = 28 and discuss the L — > 00 

extrapolation by means of the BST algorithm. We determine in a self-consistent way both 
^ I the coupling and the exponent uj of the leading correction to 74 at finite L that we assume 

■ to be ~ l/L"^. The results are 74 = -1.7422006(9) and |a; - 1| < 10"^. 
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1. Introduction 

The effort dedicated to the confirmation of Sen's conjectures [|| has established Witten's 
string field theory Q as a distinguished calculational framework where quantitative pre- 
dictions can be readily obtained. 

As is well known, Sen's conjectures concern the properties of the non perturbative 
vacuum of bosonic open string theory and relate, in the simplest case, the tachyon insta- 
bility to the presence of the D25-brane. In brief, the conjectures predict that, when the 
tachyon condensates, a non perturbative vacuum is reached whose physics is described by 
closed string theory with no tachyon and branes at all. The matching between the energy 
of the tachyon false vacuum and the D25-brane tension leads to a clear-cut prediction for 
the value of the tachyon potential at its global minimum. Therefore, the check of Sen's 
conjectures requires a reliable calculation of the potential itself. This is a typical off-shell 
quantity and string field theory is an ideal framework for the calculation. 

Unfortunately, an exact expression for the non perturbative vacuum is not known 
at present and some kind of analytic approximation is needed. In string field theory a 
possibility is offered by the level truncation approach (LT) that is a systematic expansion 
in string field components with increasing masses. The accuracy of the expansion, in a 
fixed gauge, is governed by a single parameter L related to the Virasoro level of the states 
included in the calculation. 

The application of LT to the verification of Sen's conjectures has been quite successfull 
up to date. Also, the same technique has been exploited in order to collect useful numerical 
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insights for the identification of the correct vacuum string field theory 1^] . After the initial 
and later revised calculations of the tachyon potential at level L = 4 the check of Sen's 
predictions has been extended to level L = 10 Q and more recently up to L = 18 Q. A 
good convergence with respect to L is observed and the conjectures are verified numerically 
with amazing agreement. 

Some technical questions, however, remain open. Level truncation breaks the gauge 
invariance of Witten's action and gauge fixing is not necessary. However it is usual to 
perform calculations in the Feynman-Siegel gauge that introduces a few simplifications 
(see Q for an exploration of other gauge choices). At the numerical level, the LT series in 
Feynman-Siegel gauge turned out to approach its asymptotic value monotonically at least 
in the first studies. A surprise came from the extended calculation in Q that revealed over- 
shooting at finite L in the values of the minimum of the tachyon potential. This observation 
raised some doubts concerning the best and the correct way to extrapolate the results to 
L — > oo. Several strategies have been proposed based on reasonable assumptions about 
the asymptotic properties of the LT series and it appears of some interest to investigate 
extrapolation methods in a more systematic way. 

Actually, the problem of extrapolating short expansion series is common in the analysis 
of lattice models of statistical mechanics. This is particularly true in several dimensions 
where, typically, only a few terms are known in the available perturbative expansions {e.g. 
high or low temperature). Robust algorithms have been devised to deal with such cases. 
This paper will be devoted to analyze the performance of one of these methods when 
applied to string field theory. 

We shall analyze in details a specific case that is easy to treat with level expansion 
and that permits the comparison with rather accurate alternative predictions. This is the 
quartic self-coupling 74 in the tachyon effective potential. A closed analytical expression 
exists ^] and its evaluation can be reduced to numerical quadratures that are non trivial, 
but in principle straightforward. 

The best available result is 74 = —1.742(1) Q; as we shall see later, it has been 
checked with two different level truncation schemes with good agreement at least within 
the numerical error. 

Here, we shall first improve the analytical knowledge of 74 by extending the LT series 
from the known result at level L = 20 |l^ up to L = 28. This effort involves an increase 
of the number of considered states by a factor larger than 10. Since 74 is certainly not 
the most interesting quantity to be studied, the motivation of this work is that we want a 
long LT series to analyze the various extrapolation algorithms in a clear way. Of course, a 
calculation at such level for more complicated quantities like, for instance, the full tachyon 
potential would be hardly feasible and we hope to draw conclusions holding for shorter 
series too. 

As a second step in our analysis, we shall weaken as far as possible the hypotheses 
underlying the extrapolation procedure. In particular, we shall never assume a convergence 
rate with leading correction~ 1/L, an assumption usually made in the common approaches. 
We shall apply the EST algorithm to accelerate the convergence of series with general 
leading corrections ~ where uj is an unknown exponent. We shall see that the 
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algorithm is able to determine simultaneously and self consistently both the best estimate 
for 74 and the exponent uj. 

The plan of the paper is the following: in Sec. Q we give some basic definition and 
review known facts about the quartic coupling. We also present our new results for the 
level expansion of 74. In Sec. ^ we discuss the properties of the BST algorithm and in 
Sec. 1^ we apply it to the analysis of the LT series. Appendices are devoted to technical 
details and tricks that we have found useful in the symbolic manipulation at high levels. 



2. Definition and level truncation results for 74 

Witten's open string field action is well known and reads 

where Q is the BRST operator, |^) is the ghost number 1 string field. Integration and ★ 



(2.1) 



product are well discussed in the literature and the reader can consult |11| for a review. 
The gauge invariance of S can be fixed by choosing the Feynman-Siegel gauge 60 1 ^) = 0. 
As explained by Sen, the problem of tachyon condensation can be studied by restricting 
to be a superposition with constant coefficients of open string states with ghost number 
1 belonging to a certain universal space Ti. [0] defined, in the Feynman-Siegel gauge, by 

n = Span ^ Y\ L-ki IJ b-rrij W c-m |0), with ki > 2, ruj > 2, n/ = -1 or > 1 ^ , 



i=l 



1=1 



(2.2) 



where |0) is the S'L(2,M) invariant vacuum and L„, 6„, c„ are the modes of the matter 
energy-momentum tensor and ghost fields 6, c. Due to a twist symmetry only states with 
even level L = Lq + 1 are relevant. Their number is summarized in Tab. |l[ 



level 


< 10 


12 


14 


16 


18 


20 


22 


24 


26 


28 


# states 


170 


231 


496 


1027 


2060 


4010 


7611 


14130 


25697 


45890 


total 


170 


401 


897 


1924 


3984 


7994 


15605 


29735 


55432 


101322 



Table 1: Number of states at a given level in the universal subspace H. The left part of the table 
corresponds to the most accurate existing calculation of the tachyon quartic coupling (L < 20). 



The evaluation of Witten's action can be separated in the calculation of the quadratic 
kinetic part and the cubic interaction terms. The quadratic part in Feynman-Siegel gauge 
takes the simple form 

^2 = ^ l^*Q'i' = ^{'f\coLo\^), (2.3) 

where (^'| is the BPZ conjugate of |\I'). The evaluation of 5*2 thus reduces to that of a 
long series of scalar products. These, in turn, can be computed separately in the matter 
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and ghost sectors. Selection rules can be exploited to know in advance which terms are 
necessarily vanishing. If we expand the field |^') in components 

|f) = rci|o) + ^<^^|s), \s)^ci\o), \s) en, (2.4) 

s 

(r being the tachyon field) the quadratic terms take the form 

where /i is a matrix with a well definite block structure. In particular /x connects only fields 
with the same (matter + ghost) level. 

The cubic terms, instead, mix fields with different levels. For the purpose of computing 
the quartic tachyon coupling we need to consider only the term as well as the couplings 
T'^cp. We denote the sum of such contributions by 

Ss = KgT^ + gT^Y.Ciip\ (2.6) 

i 

where k = 3^/^/2^. The evaluation of S3 can be split according to the level of 93*. The full 
action is just a cubic polynomial in the space-time independent fields. 

The (reduced) linear equations of motion of {(/?} obtained by differentiation of 5*2 + 6*3 
can be solved in terms of T. We stress that these are not the full equations of motion 
because S3 contains only a subset of the cubic couplings. Substituting the solution back 
in 5*2 + 5*3 this procedure leads to the quartic tachyon potential 

ViT) = -±T' + ngT'-^g' a' Y^{„-%c,c, (2.7) 

It is conventional to define a quartic coupling 74 independent on a' and g by writing the 
quartic term as K?a' g'^^^T^. We shall adopt this notation and present our results in terms 
of 74. 

Of course, the general construction of the tachyon potential requires to solve the full 
quadratic equations of motion of all fields with respect to T. The branch of the solution 
containing the unstable vacuum is the one with fields vanishing as T ^ 0. Expanding in 
powers of T we recover the perturbative effective potential of the tachyon. For the quartic 
term, only the leading contributions are relevant and the very simple expression in Eq. ( |2.7] ) 
can be used. This is a major simplification holding for the quartic term only. 

To evaluate Eq. ( ^.71) we need and the cubic couplings {cj}. The inversion of the 
mass matrices /U is non trivial because we are interested in an exact result. We must invert 
by symbolic manipulation and this is hard with the matrices at hand. To solve this 
technical problem we have used a specific technique based on modular arithmetic and fully 
described in Appendix |^. 

For the cubic interaction terms between states in Ti we have employed the conservation 
law technique discussed in ||T^. The general term takes the form 

j ^^^^^I- = ^(^VV" (VsIN)^^^ bV^^ (2.8) 



- 4 - 



where IV3) is Witten's three string vertex]^. For the computation of 74 we just need the 
case that we write briefly as 



(2. 



a,b,c 



The action of the lowering operators on the dual vertex (V3I can be written in terms of 
positively moded operators acting on the three strings kets. The fast derivation of these 
conservation laws at arbitrary high level is described in Appendix 

The above procedure has been implemented in a Mathematica code. As a first check 
we have reproduced the various known contributions 74^^ from fields at fixed level L for all 
L < 20 [|lO[ reported in Tab. |2| (the contributions from odd level states vanish). 



L 


74,L 


L 


74,L 


2 


34 
~27 


12 


3389771245837595 
183448998696332259 


4 


1399 


14 


78109439317035676 


6561 


5853690776582965719 


6 


18016964 


16 


298001292970739836603 


215233605 


29543127065508336986430 


8 


308423449 


18 


866910757208081799449068 




6973568802 


109838392116853446081848097 


10 


623826309310 


20 


189246300010789017377006911 


22876792454961 


29830815082559550621170156697 



Table 2: Previously known LT results for the quartic tachyon coupling 74. The symbol 74^^ denotes 
the contribution to 74 from the fields with level L. 

Then, we computed the full exact contribution of states at levels L = 22, 24, 26, 28. 
The total number of states is thus ~ 10^ to be compared with ~ 8 • 10^ when only fields 
with level L < 20 are considered. Our results are 



74,22 



74,24 



74,26 



1987896641921822231631291200 
' 381520424476945831628649898809 

83592713383361636723926342390586 



19190858871614852276752718559991509 

6244227225924504591720556698980462 
1689629965870438080888011090607948075 



0.005210, 

-0.004356, 



-0.003696, 



(2.10) 
(2.11) 
(2.12) 
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12670705380847526175723657846886589222 

74 28 = -0.003175. (2.13) 

' 3990838394187339929534246675572349035227 ^ ' 



The resulting series for 74 has 14 terms and their sum is 

22 • 26863 • 36899 • 68230523 • 365310133 • 1823562056866280202643 
383. 52 -11 -13 -17 -19 -23 



-1.70028. (2.14) 



This figure can be compared with the numerical evaluation of the available closed form 
for 74. As we mentioned in the Introduction, the best result is 74 = —1.742(1) ^. The 
discrepancy between this value and the partial sum of the LT series is 3.3% at L = 20 
and 2.4% at L = 28. It is clear that the convergence is slow and seems to be roughly 
described by a leading correction vanishing as 1/L. This conjecture is supported in this 
specific case by the analysis of a different kind of level expansion, described in Q and based 
on an explicit description in terms of a reduced set of world-sheet diagrams. It would be 
important to possess a reliable procedure to extrapolate the above table and predict the 
asymptotic L — > 00 result, possibly with some estimate of the extrapolation error. Indeed, 
with a leading error of the form 1/L even a small error of about 1% at level L = 10 would 
require an unfeasible calculation at level L = 100 to reach a 0.1 % precision. 



3. The extrapolation problem: polynomial and BST algorithms 

Given the LT series, the simplest approach to determine the asymptotic value of 74 is that 
of performing a polynomial approximation in the variable 1/L. In other words, given the 
values 74 (L) = X]^2 74,' ^ certain set of levels L G £ we choose an integer < |£| and 
solve the minimization problem 

2 

= min. (3.1) 

The resulting coefficient ao(-/V) is then chosen as an estimate of 74. In principle, better 
estimates can be obtained as N is increased up to its maximum value 

The problem with this (kind of) approach is that some assumption is implicit about the 
asymptotic convergence of LT. In particular, the expansion in powers of 1/L is reasonable 
and supported by numerical data, but could nevertheless be misleading if LT convergence 
turned out to be dominated by leading corrections of the form 1/L'^ with an exponent close 
but not equal to one. 

This effect can be seen if we can analyze a simple artificial example consisting of a 
short sequence of 10 pairs 

{(/ln,/(/ln))}n=l,2,...,10' 

{M = {l,^,^,---,^}, /(/i) = l + /i^ (3.2) 

^The non negligible error is due to the singularities affecting the numerical quadratures required for the 
calculation. 



E 



74 (L) 



N 
ra=0 



,L- 
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To determine numerically the limit /(O) = 1, we can perform a polynomial fit as in Eq. ( |3.1] ) 
with various N. For the exponent a we consider the following four values close to one: 
a = 0.8,0.9,1.1,1.2. 

In Fig. |l] we plot the estimate ao as a function of N for the four considered a. We 
see that as N increases the estimate of the asymptotic constant term improves, but at the 
largest possible A'^ there is still some non negligible disagreement. 
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Figure 1: Accuracy of the polynomial extrapolation. 

It is difficult to obtain the exact value {i.e. /(O) = 1) looking at the curves, especially 
when a < 1. Even for a = 0.9 the error at the last point is of the order of the percent. 

On general grounds, the extrapolation of the asymptotic value /oo of a sequence {/n} 
is a relatively easy problem for linearly convergent sequences defined by the condition 

fn+l foo 



\p\ < 1, with p = lim 



fn - fo 



(3.3) 



In this case rigorous theorems |13] assure the existence of algorithms able to accelerate the 
asymptotic convergence of any such sequence. Instead, when p = 1, the sequence is said to 
be logarithmically convergent and there is no general method that guarantees acceleration 
of all sequences. The best that one can do is to consider a restricted class of logarithmically 
convergent sequences with a given asymptotic form such that an accelerating algorithms 
can be found. 
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An important example of this kind, relevant to our discussion, is that of sequences 
fn ~ /oo + ain~^^ + a2n~^^ + ■ ■ ■ where < wi < u;2 < ■ ■ ■ are positive exponents. 
The convergence ratio is indeed p = 1 and there is an accelerating procedure known in 
statistical mechanics as the BST algorithm [14|. To describe it in full generality let us 



consider a function f{h) admitting the following asymptotic expansion as /i — > 

/(/i) = /(0) + ai/i'^i +02^2 +•••, 0<a;i < tJ2 < (3.4) 

and suppose that we are given a set of pairs 

{{K, fihn))}o<n<N, < hN < Hn-i < ■ ■ ■ < ho. (3.5) 

The BST algorithm attempts to determine /(O) by building a set of new sequences whose 
convergence is accelerated with respect to the initial one. To this aim, we build a grid of 
values {fn,m}o<n<N-i<m<N Setting initially 

/n -1 = 0, fn,0 = f{hn), (3.6) 

and defining the remaining values by the following iteration with respect to m: 

hn \ f ^ /ri+l,m— 1 fn,m—l 



fn,m — fn+l,m—l ~l" {fn+l,m—l /n.m— l) 



h"n+m J \ fn+l,m—l fn+l,m—2 



1 



-1 



(3.7) 

where w is a positive real number playing the role of a free parameter. 

To apply the method, we compute the BST grid and, in particular, the final value 
fo,N-i{^) as functions of lo. The best choice uj = uj* is the one that makes the last 
generated sequences (i.e. fn,m with m = — 1, — 2, . . .) most flat. There is a certain 
amount of subjective judgement in this step, but usually, the optimal uj* is identified in a 
quite narrow window. It is common to define a measure of the "non-flatness" of the last 
sequences by constructing the sum of absolute differences between a certain set of points 
in the BST grid. In particular, we can define for a certain K the border sum of differences 

K-l K 

^Kiuj) = J2 \fj,N-l-j(.'^) - fj+l,N-2-j{uj)\ + J2 \fj,N-l-j{uj) - fj-l,N-l-j(.Uj)\, (3.8) 
j=0 j=l 

and choose lo* as the value that minimize (5i<-(a;). After this, the predicted limiting value 
is simply chosen to be /o,Ar-i(c(j*). 

As a final step, one would like to determine within reasonable bounds the prediction 
error. To this aim, the quantity SxiuJ*) itself can be used and, in practice, it works rather 
well. However, we strongly suggest to inspect the grid of fn,m values and analyze the 
flatness of its subsequences in order to check this estimate. 

We now repeat the analysis of our toy data set with the BST algorithm in order to 
appreciate the advantages compared with the polynomial fit. We estimate the error using 
Eq. ( ^.^ ) with K = 6. The results are shown in Fig. ^ where we plot the difference 6e{uj) 
and the prediction /o,9(<^) as lo is varied in a wide range. For each considered a we can 
check that 5q = when lo = a. Also, as a further check of the consistency of the algorithm. 
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we can verify that 6q vanishes at co = a/2 too. The correct value u; = a is thus selected 
automatically by the algorithm. Concerning the prediction, in the second part of Fig. ^ we 
show that /o,9(a) = 1 within machine precision. These properties are actually guaranteed 
for the BST algorithm that produces rigorously a flat sequence after just two iterations 
when applied to data following the simple law Eq. ( [3.21 ). Fig. ^ just gives a numerical check 
of these facts. 




Figure 2: Accuracy of the BST extrapolation. 



The analysis of this simple example leads therefore to the following important conclu- 
sions. First, the polynomial extrapolation can be apparently converging to stable values, 
but if the leading correction is ~ 1/ L'^ with a mismatch u ^ 1 it can miss the correct result 
by a small but non zero amount. Second, the BST algorithm finds consistently both the 
exact exponent of the leading correction and the asymptotic value. It seems thus useful 
to examine the perturbative series obtained with level truncation by means of the BST 
algorithm (or variants) to confirm and possibly refine the polynomial results. 

4. Determination of 74 

After this preliminary and general considerations, we consider the LT series for 74 and 
repeat the above analysis. In order to appreciate the information associated with the 
new computed terms we first perform the polynomial fit analysis on the first 10 terms 
and separately on the full extended series. The results are summarized in Tab. ^; all the 
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shown digits are meaningful since the LT terms are known exactly and the polynomial 
minimization can be performed with no numerical approximations. 

Both sequences appear to converge well to stable values. In the first column (fit to 10 
partial sums) the difference between the last terms (A^ = 8,9) is ~ 10^^. In the second 
column (fit to 14 partial sums) this difference (now between = 12, 13) is instead reduced 
to ~ 10^''. Therefore, the extended series affords an estimate of 74 roughly one order of 
magnitude more precise in spite of the fact that the raw partial sum of the first 14 terms is 
only ~ 28/20 times closer to the predicted —1.742(1) with respect to the one with 10 terms. 
This means that the polynomial fit procedure is effective in accelerating the convergence. 

What about the BST algorithm ? Is it able to accelerate even more the convergence 
rate ? To investigate this question, we apply it to the full sequence with the conservative 
choice K = 6. The results are shown in Tab. |^ where we collect the best exponent uj and 



degree 


L < 20 


L < 28 


4 


-1.742020923 


-1.742068107 


5 


-1.742116782 


-1.742143787 


6 


-1.742172515 


-1.742181488 


7 


-1.742188837 


-1.742192207 


8 


-1.742193921 


-1.742195776 


9 


-1.742196399 


-1.742197548 


10 




-1.742198539 


11 




-1.742199101 


12 




-1.742199433 


13 




-1.742199638 



L< 




74 (BST) ±,56 


20 


0.99896 


-1.74221(6) 


22 


0.99898 


-1.74221(2) 


24 


0.99965 


-1.742202(3) 


26 


0.9998 


-1.7422009(9) 


28 


0.999875 


-1.7422006(3) 



Table 3: Results from the polynomial fit 
with data i < 20 or L < 28. 



Table 4: Results from the BST algorithm. 



the prediction for 74 when the data at level < L are used. From the table, we obtain the 
values 74 = -1.7422006(3) and |a; - 1| < 10"^. The quality of the convergence can be 
deduced by a look at Tab. |5| where we show in the 5 columns the last sequences generated 
by the algorithm in the analysis of the L = 28 series. The conclusion is that BST is able 
to achieve two important results: first, the value lj = 1 is strongly supported placing on 
firmer ground the underlying asumption made in the polynomial approach. Second, 74 is 
estimated with an even better accuracy as can be seen looking at the differences between 
the boldface figures of Tab. |5|. 

As a final comment, in order to assess the reliability of the prediction we can crosscheck 
it with another extrapolation technique as this is the safest procedure to determine the 
actual estimate error. A standard method for linearly convergent sequences and also for 
logarithmic ones with integer exponents (here a; ~ 1) is Levin u-transform |l^. We shall 
not discuss in details this algorithm, but just quote the result that is obtained in the imple- 
mentation described in ||l^. The application to the L = 28 series gives 74 = —1.7422009(6). 
The two results are completely compatible and we quote the BST central value with a very 
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iter.: 9 


10 


11 


12 


13 


1.7422001494 


1.7422003937 


1.7422005398 


1.7422005287 


1.7422005986 


1.7422003792 


1.7422004666 


1.7422005294 


1.7422005175 




1.7422004648 


1.7422004751 


1.7422005189 






1.7422004966 


1.7422005362 








1.7422005072 











Table 5: Sample subset of the BST grid. The boldface figures are would be involved in the 
calculation of S4. We show only 10 digits, but the calculation has been done with much higher 
precision. 



conservative error estimate as our final number 

74 = -1.7422006(9), |u; - 1| < 10"^. (4.1) 



5. Conclusions 



In this paper we have discussed the convergence of level truncation in open bosonic string 
field theory focusing on the specific test case of the quartic tachyon coupling. We have 
pushed the computation of this quantity up to level L = 28 and we have analyzed the 
resulting level expansion series by polynomial fits and the extrapolation BST algorithm. 
We have given numerical results supporting strongly the value uj = 1 for the exponent of the 
leading corrections at finite L of the form ~ . Also, we have obtained a very precise 
determination of the coupling itself. The proposed methodology is of general purpose and 
can be applied to any level truncation series. 
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A. Modular arithmetic and Dixon's algorithm 

In the calculation of the contributions to 74 from fields at a given level, we have to solve 
a linear problem associated to the mass matrix ^ = {^ij} to evaluate the last term in 



Eq. (2.7). It is interesting to note that /ijj G Z. Indeed, the matrix elements of are 



factorized 



M.,=/^r''^-^fr*. (A.i) 



The term //^j^"*^* can only be 0,ibl due to the ghost algebra. The term ^'^^^^^"^ is integer 
because it arises computing the scalar products of Virasoro descendents and the only dan- 
gerous (non integer) contribution can be the central charge term c/12 n(n^ — 1). However, 
if c = 26 it is evident that ^(n - l)n(n + 1) e N. 

The components of the vector of known terms in the linear problem are the cubic 
couplings Cj defined in Eq. (p.6|). Apart from a trivial \/3 factor, these are rational numbers. 
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Indeed, the factor \/3 arises in the basic 3 point ghost vertex (V3|ci|0)®ci|0)(E>ci|0) and ah 
the other manipulations do not introduce irrational numbers. Multiplying by the largest 
denominator appearing in \/3 q the problem can therefore be reduced to a totally integer 
problem Ax = b, with ^ G Z"^" and 6 G Z". 

It is known that the exact symbolic solution of such a problem can be quite hard. In 
floating point arithmetics, the complexity of solving a linear problem is only polynomial 
(cubic in the matrix dimension). However, in infinite precision it is typical to generate 
very large integers in the intermediate steps of the solution and the complexity becomes 
exponential. The problem is not trivial at the level we are working where the dimension of 
the relevant matrices can be as large as ~ 700 with integer entries with several tenths of 
digits. 



A simple solution to this computational problem is afforded by Dixon's algorithm |18|, 
a procedure arising in the field of modular p-adic arithmetic. Suppose we want to find the 
rational solution of 

Ax = b, ^gZ"^'', 6gZ". (A.2) 

We solve the problem in 3 steps. In the first step we choose a prime p such that the 
following matrix exists: 

C = A^^ mod p. (A.3) 

The calculation of C is not computationally expensive because all integers manipulated in 
the inversion of A are bounded by p. In the second step we choose an integer m and define 

6(0) = 6; = Cb^^ mod p, b'^'^^^ = p-\b^^ - Ax^''>), < i < p. (A.4) 

At the end of the iteration we set 

m— 1 

x=Y^ pVW G Z^. (A.5) 

2 = 

It can be shown that x is a solution of the initial problem mod p'^. Finally, in the third 
step, we reconstruct the rational solution by solving the p-adic representation of x. For 
each 1 < J < n we define 

^(-1) = ^(0) ^ ^.^ ^(-1) ^ 0, = 1. (A.6) 

Denoting by fioor(a;) the greatest integer less than or equal to x, we start with i = and 
iterate 

(A.7) 



while(n(*) 


> p"/2){ 


Q 


= fioor(u(^"^V^^*^) 










i 


i + l, 
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At the end, we assign the j-th component of the proposed solution according to 

x,{p,m) = {-iy-^ (A.8) 

where i takes the last value attained during the cycle in Eq. ( |A.7| ). The above procedure 
depends on the initial prime p and the integer m as emphasized in the notation. It can 
be shown that for any p, x{p, m) equals the unique rational solution of the linear problem 
problem as soon as m exceeds a value that is in principle unknown, but that is assured 
to exist. Therefore, in practice, one chooses a reasonably large p (the larger p the more 
expensive is the determination of C) and than increases m until x(p, m) satisfies the initial 
system, a check that is not computationally expensive. 



B. Conservation laws at high levels 

As explained in |12|, the evaluation of the cubic vertex among states belonging to Ti can 
be reduced to algorithmic rules that are fixed in terms of a set of meromorphic vector fields 
and quadratic differentials satisfying suitable conditions. Without repeating the details 
of the construction, we just give a few simple formulas that permit a totally automatic 
determination of these objects with a negligible computational effort. We hope that this 
Appendix will be useful for the reader interested in reproducing or extending this kind of 
calculation. 

The discussion can be done separately for the matter and c{z) ghost sectors. The 
conservation laws for the h{z) ghost are closely related to those in the matter sector. We 
define 

and expand in powers of z the following maps 

, , , /2 \ 2 10 o 38 . 574 7 30050 q , 

f2{z) = tan - arctan z] = -z z'^ -\ z^ z' -\ z^ ^ , (B.2) 

^ ' V3 y 3 81 729 19683 1594323 ' 

f^{z) = S{h) = -h{-z). (B.4) 
B.l Matter sector 

Following the notation of [|l^], the conservation laws for the Virasoro operators are deter- 
mined by vector fields Vn{z), n E N of the form 

Vn{z) = ^j^K(^), (B.5) 

where Vn is a. polynomial with degree [(n — l)/2] {[x] denotes the integer part of x) and 
the following expansions around z = hold 

Vn{f2) 1 , / 0{z) n even VnUifi) ,^ „^ 

= + 0{z-) n odd ' = (^-^^ 
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In 1 12 1, the fields f„ are built in sequence starting from vi and defining each of them in terms 
of suitable linear combination of the previous fields. This requires to fix the coefficients of 
the linear combination at each step. The whole procedure is not very expensive, but we 
think that it is interesting to give a compact formula for the result. Solving to the desired 
order the constraints obeyed by f„ we find: 



where F{z) = tan(3/2 arctan z) and Kf{z) is defined in terms of the Laurent expansion of 

M 
f{z) = CnZ^ ^ K{f{z)) = CnZ^. (B.8) 

n=-N n=-N 

The application of the algorithm determines very quickly all the vector fields Vn- To give 
an example, let us consider n = 10. We have 

1 1 



z2-3F'(z)F(z)9 



(B. 



_ -1024 14336 _ 99136 461120 _ 1668820 ^ 

~ 177147 z9 ^ 531441^7 ~ 1594323 ^ 4782969^3 ~ 14348907 z ^ 

and thus acting with K and multiplying by 2;2 — 3 we find 

4 (20736 - 96768 + 223056 - 345840 z^ + 417205 z^) 

Vioiz) = i '-. (B.IO 

' 14348907 ^ ' 

The list of the first 9 vector fields is 

IV ; g, y ■^^ 81 

8 (-6 + 13^2) 2 (48- 124^2 + 163 z^) 
729 ' = 5187 

4 (48- 144 z2 + 217 z^) 2 (-576 + 1968 z^ - 3352 z'' + 3967 z^) 

^ ^ 6561 ^ ' 59049 

16 (-144 + 552 z2 - 1050 z"^ + 1367 z6) 



Vsiz) = 



2 (2304 - 9792 z2 + 20592 z^ - 29336 z^ + 32969 z^) 



531441 
B.2 Ghost sector 

In the ghost sector, the conservation laws for the modes of the c(z) ghost field require the 
introduction of quadratic differentials (pn{z){dz)'^ with general form 
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where ^n{z) is a polynomial of degree 2[n/2]. The additional conditions that fix are 



Mf2)if'. 



1^2 1 I j ^(1) even 
0{l/z) n odd 



+ 



<Pn{h,z){f[.Y = 0{l/z). 



{B.m 



To find an explicit solution we compute 

F'{zf 



K[{z'-?>) 



Cn+2 Cn+1 _i_ £|i 



where p = 0(1) if n is even(odd). The solution is then simply 



(z2-3)</.„(z) 



Cn+2 Cra+1 
^n+2 z^~^^ 



+ ■■■ + 



Again, we give an example. For n = 4 we have 



, 2 N F"(z) 16 16 26 200 n 



(B.13) 



(B.14) 



(B.15) 



'F(z)6 27 z6 81 z4 81 z2 729 
and we obtain (/>4 dividing the sum of the first 3 terms by z'^ — 3. The first 10 cases are 

-2 (2 + z2) ^ ^ _2 (l2 + z2) 



$o(^) = -3, <^i{z) 



, , 2 (-24 + 8^2 + 13^4) 



3 ' ^'^""^ 27 
2 (-48 + 36 22 + 172"^) 



$6(^) 



$7(z) 



81 ' ' 243 

-2 (288 - 336 ^2 + 2 z'' + 217 z^) 
2187 

-2 (576 - 912 z2 + 312 + 391 z^) 
6561 



^9(2) 



2 (-1152 + 2304^2 _ 1440^4 _ 416^6 + 1367 2^) 
19683 

2 (-20736 + 50112 z'^ - 44208 z'^ + 5640 z^ + 25703 z^ 

531441 



«'io(^) = 
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